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Abstract. We present optical and near-infrared observations with Keck of the binary millisecond pulsar PSR 10751 + 1807. We 
detect a faint, red object - with R = 25.08 ± 0.07, B - R = 2.5 ± 0.3, and R - I = 0.90 ± 0. 10 - at the celestial position of the 
pulsar and argue that it is the white dwarf companion of the pulsar. The colours are the reddest among all known white dwarfs, 
and indicate a very low temperature, T c ff ss 4000 K. This implies that the white dwarf cannot have the relatively thick hydrogen 
envelope that is expected on evolutionary grounds. Our observations pose two puzzles. First, while the atmosphere was expected 
to be pure hydrogen, the colours are inconsistent with this composition. Second, given the low temperature, irradiation by the 
pulsar should be important, but we see no evidence for it. We discuss possible solutions to these puzzles. 
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1. Introduction 

Among the pulsars in binaries, the largest group, the low-mass 
binary pulsars, has low-mass white-dwarf companions. Before 
the companions became white dwarfs, their progenitors filled 
their Roche lobe and mass was transferred to the neutron stars, 
thereby spinning them up and decreasing their magnetic fields. 
Considerations of the end of this stage, where the white dwarf 
progenitor's envelope becomes too tenuous to be supported 
further, allow one to make predictions for relations between 
the orbital period and white dwar f mass, and orbital period 
and eccentricity (for a review, e.g., Phinnev & Kulkarni 1994; 
IStairsll2004l) . Furthermore, after the cessation of mass trans- 
fer, two clocks will start ticking at the same time: the neutron 
star, now visible as a millisecond pulsar, will spin down, while 
the secondary will contract to a white dwarf and start to cool. 
Consequently, the spin-down age of the pulsar should equal the 
cooling age of the white dwarf. 

From optical observations of white-dwarf companions 
to millisecond pulsars one can estimate the white-dwarf 
cooling age and compare it wi th the pulsar spin-down age . 
Initial attemp ts to do this ( lHansen & Phinnev 1998a b; 
ISchonberner et al.ll200C4) revealed a dichotomy in the cooling 
properties of white dwarfs in the sense that some white 
dwarf companions to older pulsars have cooled less than 
those of younger p ulsars. In particular, the companions of 
PSR J0437-4715 JPanzieer et alJ 1 19931 Ivan Straten et all 
2001b and PSR B 1855+09 Jvan Kerkwiik et alF koOO ; 
Rvb a & Tavlorl Il99ll) have temperatures of about 4000- 
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5000 K, with characteristic pulsar ages of 5 G yr. This is in 
contrast to the companion of PSR J1012+5307 JLorimer et alJ 
1199.4 Ivan Kerkwiik et alJll99fil ICallanan et alJll998l) . which 
has a higher temperature (8600 K), while it orbits an older 
pulsar (8.9 Gyr). 

A likely cause for this dichotomy is the difference in the 
thickness of the envelope of hydro gen surroundi ng the he- 
lium core of the white dwarf jAlberts et alJll996l) . After the 
cessation of mass transfer, the white dwarfs have relatively 
thick (~ 1O~ 2 M ) hydrogen envelopes which are able to sus- 
tain residual hydrogen shell-burning, keeping the white dwarf 
hot and thereby slowing the cooling ( lDriebe etalJll998l) . The 
shell burning, however, can become unstable and lead to ther- 
mal flashes which can reduce the mass of the envelope. White 
dwarfs with such reduced, relatively thin (5; 1O~ 3 M ) hydro- 
gen envelopes cannot burn hydrogen and, as a result, cool 
faster. The transition between thick and thin hydrogen en- 
velopes was predicted to lie near 0.1 8-0.20 M R (where heav- 
ier white dwarfs have thin envelop es; lAlberts et alJ [l996; 
ISarna et allipOollAlfhaus et all200ll) . 

Until recently, PSR J1012+5307, with an orbital period 
Pb - 0.60 d, was the only system for which a thick hydro- 
gen envelope was required to match the two timescales. Given 
the relati on between the white dwarf mass and the orbital pe- 
riod (|Joss et al.ll987t iRappaport et alJl995t iTauris & Sav onii e 
1999), companions in similar or closer orbits should have 
similar or lower mass, and thus have thick hydrogen en- 
velopes as well. This was confirmed by the recent discovery 
of two new, nearby, binary millisecond pulsars with orbital pe- 
riods similar to that of PSR J1012+5307; PSR J1909-3744 



2 



C. G. Bassa et al.: The ultra-cool white dwarf companion of PSR J075 1 + 1807 



(1.53d, I.TacobvetalJl2Q05h and PSR J1738+0333 (0.354d, 
Jacoby et al., in prep.; see Ivan Kerkwiik et ai1l2005l for pre- 
liminary results). For both, the temperatures and characteris- 
tic ages are similar to those of PSR J1012+5307, and thus 
one is led to the same need for a thick hydrogen envelope. 
These discoveries, combined with the thin envelopes inferred 
for PSR J0034-0534 ( 1 .59 d) and binaries with longer periods, 
suggest that the transition occur s at a mass that correspon ds to 
an orbital period just over 1.5 d ( Ivan Kerkwiik et all2005t) . All 
systems with shorter orbital periods should have thick hydro- 
gen envelopes. 

The two known millisecond pulsars with white dwarf com- 
panions that have shorter orbital periods than PSR J1012+5307 
but do not have optical counterparts, are PSR J0613-0200, 
with a 1.20d period, and PSR J075 1 + 1807, which has the 
shortest orbital period of all binary milli second pulsars with 
M c > 0. 1 M companions, 0.26 d ilLundgren et a l J l 199.4 . The 
latter system is of particular interest because the companion 
mass has been determined from pulse timing (Mwd = 0.19 + 
0.03 M at 95% confidence: lNice et all2005l) so that one does 
not have to rely on the theoretical period-mass relationship. 
Intriguin gly, for PSR J075 1 + 1 807, optical observations from 
iLundgren et alJ ll 19961) set a limit to the temperature of 9000 K, 
which is only marginally consist ent with it having a thick hy- 
drogen envelope. Based on this. lErgma et alJd200ll) . suggested 
the hydrogen envelope may have been partially lost due to irra- 
diation by the pulsar. 

The faintness of the companion to PSR J075 1 + 1807 
aroused our curiosity and motivated us to obtain deep observa- 
tions to test the theoretical ideas discussed above. We describe 
our observations in Sect. [5] and use these to determine the tem- 
perature, radius and cooling history in Sect. [3] In Sect. |4] we 
investigate irradiation by the pulsar, finding a surprising lack 
of evidence for any heating. We discuss our results in Sect. [5] 

2. Observations and data reduction 

The PSR J075 1 + 1807 field was observed with the 10 me- 
ter Keck I and II telescopes on Hawaii on five occasions. On 
December 11, 1996 the Low Resolution Imaging Spectrometer 
(LRIS, lOke et alJ 1 19951) was used to obtain B and fl-band 
images, while the Echellette Spectrograph and Imager (ESI, 
ISheinis et alJl2002T) was used on December 2 1 , 2003 to obtain 
deeper B and R -band, as well as /-band images. The /?-band fil- 
ter used that night was the non-standard "Ellis R" filter. The ob- 
serving conditions during the 1996 night were mediocre, with 
078-171 seeing and some cirrus appearing at the end of the 
night. The conditions were photometric during the 2003 night, 
and the seeing was good, 076-078. The third and fourth visit 
were with LRIS again, now at Keck I, on January 7 and 8, 
2005. The red arm of the detector was used to obtain /?-band 
images. The seeing on the first night in 2005 was rather bad, 
about l'/5 and improved to about 170 on the second night. The 
conditions on these nights were not photometric. Finally, a se- 
ries of 36 dithered exposures, each consisting of 5 co-added 
10 s integrations, were taken through the filter with the Near 
Infrared Camera (NIRC: lMatthews & Soiferlll994 on January 
26, 2005. The conditions were photometric with 076 seeing. 



Table 1. Observation log. 



Field Time (UT) 


Filter 


Mnt V a / 


sec z 


December 11, 1996, LRIS 










r 


2+10 


1.07 


08-27-08-29 


g 


2+10 


1.07 


SA9*S 09-78-09-31 


g 


2+10 


1.07 


09-33-09-3*S 
u™. J ~j \jy . j j 




2+10 


1.08 


PSR J075 1 + 1807 09:45 


R 


10 


1.39 


09:47-09:59 


R 


2x300 


1.36 


10:01 


R 


600 


1.31 


10:13 


B 


600 


1.26 


December 21, 2003, ESI 








PSRJ0751 + 1807 10:06-10:27 


R 


3x360 


1.14 


10:29-10:57 


I 


6 x 240 


1.08 


11:00-11:33 


B 


3 x 600 


1.04 


NOC2419 11:40 


B 


10 + 30 


1.06 


11:44 


R 


10 + 30 


1.06 


11:47 


I 


10 + 30 


1.06 


January 7, 2005, LRIS 








PSRJ0751 + 1807 11:54-12:53 


R 


5 x600 


1.05 


January 8, 2005, LRIS 








PSRJ0751 + 1807 11:42-12:51 


R 


6x600 


1.05 


January 26, 2005, NIRC 








PSR J075 1 + 1807 08:06-08:56 


A' s 


36x50 


1.07 


2MASS star 3 08:59 




0.4 


1.02 



2M AS S J075 1 062 1 + 1 807253 



Standard stars llLandoltlll99i IStetsonl [20001 were observed 
in 1996 and 2003, while a 2MASS star JCutri etai1l2003l) in 
the vicinity of PSR J075 1 + 1807 was observed to calibrate the 
NIRC data. A log of the observations is given in Tabled 

The images were reduced using the Munich Image Data 
Analysis System (MIDAS). The BRI images were bias- 
subtracted and flat-fielded using dome flats. The longer expo- 
sures in each filter were aligned using integer pixel offsets, and 
co-added to create average images. The near-infrared images 
were corrected for dark current using dark frames with identi- 
cal exposure times and number of co-adds as those used for the 
science frames. Next, a flatfield frame was created by median 
combining the science frames. After division by this flatfield, 
the science frames were registered using integer pixel offsets 
and averaged. 

2.1. Astrometry 

For the astrometric calibration, we selected 14 stars from 
the secon d version of the USNO CCD Astrograph catalogue 
(UCAC2: Zacharias et alJl2004 that overlapped with the 10 s 
/?-band LRIS image of December 1996. Of these, 1 1 were not 
saturated and appeared stellar and unblended. The centroids 
of these objects were measured and corrected for geometric 
distortion using the bi-cubic function determined by J. Cohen 
(1997, priv. comm.) 1 . We fitted for zero-point position, plate 
scale and position angle. The inferred uncertainty in the single- 

1 http://alamoana.keck.hawaii.edu/inst/lris/coordinates.html 
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star measurement of these 1 1 stars is 07057 and (X'083 in right 
ascension and declination, respectively, and is consistent with 
expectations for the UCAC measurements of approximately 
07020 for stars of 14th magnitude and (X'070 for stars 2 magni- 
tudes fainter. 

This solution was transferred to the lOmin 7?-band LRIS 
image using 91 stars that were present on both images and were 
stellar, unsaturated and not blended. Again the zero-point posi- 
tion, plate scale and position angle were left free in the fit and 
the final residuals were 07016 and 07019 in right ascension 
and declination. The UCAC is on the International Celestial 
Reference System (ICRS) to £ (X'01, and hence the final sys- 
tematic uncertainty with which our coordinates are on the ICRS 
is dominated by our first step, and is ~ 0703 in each coordinate. 

Our images, with the position of PSR J0751 + 1807 
jNice et al]l2005ll indicated, are shown in Figure ^ On the 
lOmin LRIS tf-band images from 1996 and 2005, we find 
a faint object, hereafter star X, at the position of the pul- 
sar. It is also, though marginally, present in the two 5 min R- 
band images from 1996, but not detected in the 10 min 73-band 
LRIS image of that observing run. Star X is clearly present in 
the 2003 ESI R and /-band images, and marginally in the B- 
band image. It is not detected in the near-infrared observations 

(Fig.nj. 

Positions for star X and other objects inferred using the as- 
trometry of the 10 min LRIS 7?-band image are listed in Table|2] 
The pulsar po sition at the time of the 1996 LRIS observa- 
tion, using the iNice et al J (|2005) position and proper motion, 
is Qrjaooo = 07 h 51 m 09 s . 1574(1), 6 ]2 (m = +18°07'38'.'624(10). 
We find that star X is offset from the pulsar position by 
-0701 ± 0706 in right ascension and 0704 + 0706 in declina- 
tion, well within the lcr uncertainties (including those on the 
pulsar position). Given the low density of about 47 stars per 
square arcminute and the excellent astrometry, the probabil- 
ity of a chance coincidence in the 95% confidence error cir- 
cle, which has a radius of 0724, is only 0.1-0.2%. Since, as 
we will see, it is hard to envisage how the companion could be 
fainter than the object detected, we are confident that star X is 
the companion of PSR J075 1 + 1 807. 

2.2. Photometry 

The DAOPHOT II package JStetsonl fl987l) . running inside 
MIDAS, was used for the photometry on the av eraged images. 
We followed the recommendations of IStetsonl lll987ll : instru- 
mental magnitudes were obtained through point spread func- 
tion (PSF) fitting and aperture photometry on brighter stars was 
used to determine aperture corrections. 

For the calibration of the optical images, instrumental mag- 
nitudes of the standard stars, determined using aperture pho - 
tometry, were compared against the values of IStetsonl (12000). 
We used the standard Keck extinction coefficients of 0. 1 7, 0. 1 1 
and 0.07 mag per airmass for B, R and /, respectively. Colour 
terms were not required for the LRIS B and R bands, but were 
significant for the ESI bands: 0.107(5 - R) for B, 0.083(73 - R) 
for R, and -0.004(7? - 7) for 7, i.e., the ESI B, R are redder 
than the standard bands, while ESI 7 is slightly bluer. The 
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Fig. 1. Images of the field of PSR J075 1 + 1807. The upper fig- 
ure shows a 1' x 1' subsection of the averaged 6 x 10 m 7?-band 
image obtained with LRIS on January 8, 2005. The bottom four 
figures show 20" x 20" subsections of the B, R and /-band av- 
erages observed with ESI in 2003 and the 7f s -band image ob- 
served with NIRC in 2005. The tick marks all have a length of 
1" and are centered on the pulsar timing position. 



root-mean-square residuals of the ESI calibrations are about 
0.05 mag in B, and 0.03 mag in R and /, while those of the LRIS 
calibration are 0.08 mag in B and 0.05 mag in R; we adopt these 
as the uncertainty in the zero-points. The near-infrared obser- 
vations were calibrated through aperture photometry with 175 
(lOpix) apertures using the 2MASS star, fitting for a zero-point 
only, as the difference in airmass between the science and cal- 
ibration images is small. We adopt an uncertainty in the Ks 
zero-point of 0. 1 mag. 

Calibrated ESI magnitudes for star X and selected other 
stars in the field are listed in Table |2 Star X is barely above 
the detection limit of the ESI 73-band observations, hence the 
large error. It is not detected in the LRIS fi-band and the NIRC 
7T s -band observations, and, scaling from the magnitude of a star 
with a signal-to-noise ratio of about 10 and 6, we estimate the 
3cr detection limits at B — 26.8 and 7f s = 21.3, respectively. 
The former is consistent with the ESI detection. None of the 
stars in Tableware covered by the small 38" x 38" field-of- 
view of NIRC, hence we do not have near-infrared magnitudes 
for these. 

The 1 996 LRIS fl-band magnitude is 25. 13 ± 0.1 1 , which is 
consistent with the ESI measurement. Since the conditions dur- 
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Table 2. LRIS Astrometry and ESI photometry of the compan- 
ion of PSR J0751 + 1807 and stars in the field. The nomencla- 
ture of the stars is according to Fig.^ The uncertainties listed 
in parentheses are instrumental, i.e., they do not include the 
zero-point uncertainty in the astrometric tie (about 0703 in each 
coordinate) or of photometric calibration (0.05 mag in B and 
0.03 mag in both R and /). 



ID 


(22000 
h m s 




#2000 B 

O / ft 


R I 


X 


07 51 09.158(4) 


li 


! 07 38.66(6) 27.56(25) 25.08(7) 24.18(7) 


A 


07 51 09.933(1) 


li 


107 05.97(1) 21.73(1) 


19.30(1) 18.31(1) 


B 


07 51 10.844(1) 


If 


107 52.91(1) 22.80(1) 


21.03(1) 20.32(1) 


C 


07 51 10.891(1) 




1 07 35.69(1) 24.30(2) 


21.81(1) 20.63(1) 


D 


07 51 10.739(1) 




1 07 32.79(1) 24.28(6) 


22.50(5) 21.99(6) 


E 


07 51 08.519(1) 


U 


1 07 59.89(2) 24.56(7) 


22.87(5) 22.38(8) 


F 


07 51 08.859(2) 


li 


107 08.83(3) 24.94(4) 


24.00(5) 23.29(4) 


G 


07 51 08.908(4) 


If 


1 07 35.71(5) 25.65(8) 


24.51(5) 23.85(6) 


H 


07 51 10.691(3) 


18 07 24.69(6) 25.69(7) 


24.94(9) 24.34(8) 



ing the 1996 LRIS observations may not have been photomet- 
ric, however, this may be a coincidence. To check for variabil- 
ity, we tied the instrumental LRIS R band magnitudes directly 
to the ESI R and / ones, using 38 stars that both images had in 
common and that had magnitude uncertainties below 0.1 mag. 
As expected given the non-standard "Ellis R" filter on ESI, we 
required a large colour term, -0.302(/?; nst - /i nst ), but with this 
the fit was adequate, with root-mean-square residuals of 0.14 
mag. Compared to the fit, the ESI minus LRIS difference in R- 
band magnitude is insignificant, -0.03 + 0.1 3 mag. Similarly, 
comparing instrumental /?-band magnitudes from 2005 January 
7 with those taken 2005 January 8 and 1996 December 11, 
fitting for an offset only, results in magnitude differences of 
0.03 + 0.07 and -0.16 + 0.12mag, respectively. Thus, no large 
variations in brightness are seen; we will see in Sect.0]that this 
is somewhat surprising. 

3. Temperature, radius, and cooling history 

We use our observations of star X, the companion of 
PSR J075 1 + 1807, to constrain its temperature, radius, and at- 
mospheric constituents, and discuss our result that the white 
dwarf does not have the expected thick hydrogen envelope. 

3.1. Colours, temperature, and atmospheric 
composition 

We first use the colours of star X to constrain its temperature. 
The red colours are largely intrinsic, as the maximum redden- 
ing towards PSR J07 5 1 + 1807 (I = 202.73 b = 21 .09) is small, 
E B - V = 0.05 ± 0.01 dSchlegeletalJll998l) . This value is con- 
sistent with the low value found for the interstellar absorption 
N H ~ 4 x 10 20 cnr 2 , as estim ated from ROSAT X -ray obser- 
vations of PSR J075 1 + 1 807 bvlBecker et all dl996l) . For com- 
parison, the relation bv |Prederil&^chrrirtTfl995Tr predicts an 
Ah ~ 3 x 10 20 cm' 2 for t he above reddening. Given the distance 
of ~ 0.6 kpc iNice et al]l2005l) . we expect most of the redden- 



ing to be in the foreground to the pulsar. Hence, the dereddened 
colours are (B - R) = 2.40 + 0.27 and (R - 7) = 0.86 + 0.10. 

In Fig.|5^, we compare the intrinsic colours of star X with 
those of other white-dwarf companions of millisecond pulsars, 
other white dwarfs, and models. We find that the colours of star 
X are the reddest for any known millisecond pulsar companion 
or white dwarf. The pulsar companion that comes closest is that 
ofPSR J0437-4715(ff - R = 2.12 + 0.06, R - / = 0.56 + 0.02 
llDanziger et al.l fT993ll and negligible extinction 2 ); the most 
similar white dwarf i s WD 0346+246 (B - R = 2.2 + 0.1, 
R- I = 0.76 + 0.08, Qppe nheimer et alJl200lh . Thus, star X 
is likely as cool or even cooler than the T W 3 700 K inferred 
for those two sources (PSR .T0437-4715: lDanziger et alJl993l: 
Hansen 2002, priv. co mm.; WD 0346+246: Qppenhe imer et alJ 
l200lUBergeronl200ll) . 

Also shown in Fig. Ek are colour s expected from model 
atmospheres of ISerenelliet al] (12001 h and of Hansen (2004, 
priv. comm.), which are specifically tailored to the low-mass, 
helium-core companions of millisecond pulsars, as well as 
those for updated low- gravity (logg = 7), pu re hydrogen at- 
mosphere models 3 of Ber geron et al.l i 19951) . One sees that 
the colours of the companion of PSR J0437-4715, as well 
as those of the hotter companions of PSR J1012+5307 and 
J0218+4232, are consistent with these models. For star X, how- 
ever, the colours are not consistent, as the models never venture 
redwards of R - I * 0.7 and B - R * 2.0. 

The change in direction of the tracks is seen in all models 
for hydrogen-rich, metal-free atmospheres; it reflects a change 
in the dominant source of opacity, from bound-free absorption 
of H at higher temperatures to collision-indu ced absorption 
of H? at lowe r ones dLenzuni et al.lll99H ISaumon et al.lll 994t 
lHansenll9 98^. The latter process is highly non-grey, and leads 
to absorption predominantly longward of the /?-band. As a re- 
sult, the R — I colour becomes bluer with decreasing tempera- 
tures, while B - R remains roughly constant. 

Could star X have a different composition? Due to the high 
gravity of white dwarfs, metals settle out of the atmosphere. 
However, some white dwarfs have atmospheres dominated not 
by hydrogen, but by helium. For the latter, the opacity sources 
are all fairly grey, and hence the colours continue to redden 
with decreasing temperatures. Indeed, the colours of star X are 
consistent with the predicti ons of th e updated log g — 7 pure 
helium models after Bergeron et al.l 111995) at T e g 4200 K 
(Fig. 00- 

From an evolutionary perspective, however, a pure he- 
lium atmosphere is not expected. Low-mass white dwarfs 
such as the companions to millisecond pulsars are all formed 
from low-mass stars whose evolution was truncated by mass 
transfer well before helium ignition (for recent models, see 
iTauris & Savoniielfl999llNelson et alJl2004 . As a result, they 
should have helium cores surrounded by rela tively thick, 
0.01 to 1% of the mass, hydrogen envelopes JPriebe et alJ 

2 As inferred from the dust maps of Schleg eTet alJ (1998); 
iDanziger et alJ fT993l) estimate E B -v = 0.07 from the work of[Knude 
J 19791) . 

3 For updated versions of the iBergeron et al 1 ll995t) models, see 
http://www.astro.umontreal.ca/~bergeron/CoolingModels/ 
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Fig. 2. a Colour-colour and b colour-magnitude diagram for 
the companion of PSR J075 1 + 1807, other millisecond pulsar 
companions, field white dwarfs, and model predictions. Show n 
with error bars are PSR J0437-4715 JPanziger et alJll99l . 
PSR J1012+5307 JLorimeretai]ll995l) . PSR J0218+4232 
jBassa et alJl2003h and PSR J075 1+1807 (this w ork), as well 
as the ultra-cool field w hite dwarfs LHS 3250 (Harris et all 
1 19991) WD 0346+ 246 Jopnenheim er et all l200lb . and GD 
392B JFarihil2 004). In the colour-colour diagram, also the full 
sample of field white dwarfs of Berger on et alJj200ll) is shown, 
with filled and open circles indicating white dwarfs with and 
without Ha in their spectrum, respectively. For the colour- 
magnitude diagram, we used parallax distances where avail- 
able and distances inferred fr om the pulsar d ispersion mea- 
sure otherwise; we omitted the B ergeron et alJ sample because 
of the large scatter (even though it likely is largely intrinsic). 
In both panels, the continuous light gr ey lines depict helium - 
core white dwarf cooling models by ISerenelli et all l200lh . 
with masses as indicated. The continuous dark-grey track is 
for a .2M H model by B. Hansen (priv. comm.; see descrip- 
tion in Hansen & Phinnev 1998a). The dark-grey dashed and 
dashed-dotted li nes ar e updates of the white dwarf models by 
Ber geron et alJ ll 1995b . for DA (hydrogen rich) and DB (he- 



ll998HAlthaus et alJ200lh . In deed, among the low-mas s white- 
dwarf companions to pulsars (van Kerkwiik et al .120051) as well 
as am ong low-mass white dwarfs in general JBergeron et alJ 
2001), only hydrogen-dominated atmospheres have been ob- 
served. 

In principle, at low temperatures, the hydrogen envelope 
might become mixed in with the helium core. Even if fully 
mixed, however, the remaining amounts of hydrogen would 
strongly influence the spectrum. Indeed, the effects of collision- 
induced absorption increa se with incre asing helium abundance 
up to N(Ue)/N(U) 10 5 Bergeron & Leggettl2002l) . 

From Fig.|2t, it is clear that the predictions for hydrogen- 
dominated atmospheres are also a somewhat poor match to 
the colours of the cooler normal white dwarfs with hydro- 
gen in their at mospheres (as inferred from absorption at Ha, 
Ber geron et all2001t filled circles in the figure). For most, this 
app ears to be due to missing blue opacity in the models (see 
iBergeron et aljEoOll for a detailed study); the visual through 
infrared fluxes are reproduced well by the models, and show 
unambiguously that collision-induced absorption by H2 is im- 
portant. Indeed, the absorption is evident in the optical colours 
of some objects, in particular LHS 3250 (shown in F ig.Et and 
SDSS J133739.40+000142.8 iBergeron & Legged I2OO2I and 
references therein). 

For our purposes, however, the case of the ultra-cool white 
dwarf WD 0346+246 is most relevant. For this source, the 
colours cannot be reproduced with either pure hydrogen or 
helium, but require a mixed atmosphere, dominated by he- 
lium (with fractional hydrogen abundances ranging from 10~ 9 
to 10 , depending on assumptions about the contribution 
of oth er opacity sources; Opnen heimer et all 1200 it iBergeror] 
2001, though recent work puts these abundances in to doubt, P. 
Bergeron 2005, priv. comm.). For all cases, the temperature is 
around 3700 K. The similarity in the colours of WD 0346+246 
and star X would suggest that star X has a similar, maybe 
slightly lower, temperature. 

From the above, we find that we cannot determine the 
temperature of the companion of PSR J075 1 + 1807 with cer- 
tainty, since we do not know its composition. Most likely, 
however, it is somewhere between the temperature inferred for 
WD 0346+246 and that indicated by the (pure helium) models, 
i.e. in the range of, say 3500-4300 K. 

A more stringent test could be provided by the near-infrared 
observations, as the R - K colour (which is similar to R - K s ) 
differs for different predictions. At a temperature of 4000 K the 
log g — 7 Ber geron et al.l Jl995l) models predict R - K colours 
of 2.7 and 1 .6 for pure helium and pure hydrogen atmospheres, 
respectively. For the same temper ature,/?-^ =1. 6 is predicted 
by the O.196M model bv lSerenelli et al.1 d200ll) . Finally, for 
WD 0346 +246, with presumably a m ixed hydrogen/helium at- 
mosphere, Oppenheimer et al. (2001) observed/? — K = -0.7. 
Unfortunately, our near-infrared observations only limit the 
colour to R - K < 3.8, which does not constrain any of these 
predictions. 
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3.2. Brightness, distance and radius 

So far, we have only discussed the colours and temperature. 
We now turn to the absolute magnitude and radius. In Fig. EJ), 
we show Mr as a function of R - /. For star X, we com- 
puted the absolute 7? -band magnitude Mr using the parallax of 
n = 1 .6 ± 0.8 mas as measured through radio timing dNice et all 
The resulting distance of 0.6^2 kpc is consistent with 
that estimated from the dispersion measure which predicts 1 . 1 ± 
0.2 kpc, using a d ispersion measure of 30.2489 + 0.003 pc cm -3 
JNice et alJl2005ll and the recent model of the Galactic elec- 
tron distribution of lCordes & Lazio1d2002l) . Correcting for the 
reddening, this implies Mr - 15.97+° ^. 

Given the similarities in the above absolute magnitude 
of star X and that of WD 0346+246 (Mg = 16.1 + 0.3; 
lHamblvet al1ll999t lOppenheimer et all EoOl ), and assuming 
similar temperature, one finds that the radius of star X should 
be comparable to the R = 0.010 R for WD 0346+246 
( Ber geronl200ll) . However, the large uncertainty in the parallax 
of PSR J075 1 + 1807 allows radii between 0.007-0.021 R . For 
the white-dwarf m ass of ~ 0.19 M inferred from pulse timing 
dNice et alJkOoH) . this is consistent the ~ 0.022 R expected 
from the 0.196 M model by dSerenelli et al.l200ll) . 

As can be seen in Fig. El the absolute magnitude is also 
consistent with t he predicted values fr om the log g — 7 pure 
helium model by Berg eron et alJ d 19951) . At a temperature of 
Jeff = 4250 K, this model has a radius of 0.020 R and a mass 
of 0.15 M , somewhat smaller than the observed O.19M . To 
correct for the small difference in mass, we computed white 
dwarf radii for the observed temperature and mass of the com- 
panion and used these to scale the absolute magnitudes of the 
pur e helium track in F ig. El At O.19M and T eff = 4000 K, 
the IPanei et all (EoOO) helium core white dwarf mass-radius 
relation predicts 0.021 R . This is ver y similar to the radius 
predicted by the Bergeron et al. ( 1995) \ogg = 7 pure helium 
models, and as such, the absolute magnitudes are comparable. 
We conclude that, with in the large uncertainties on the paral- 
lax distance, the absolute magnitude and radius that we derive 
for the companion of PSR J075 1 + 1807 are consistent with the 
predictions for a pure helium atmosphere. 

We note that of t he models presented in Fig. E] those of 
Ber geron et al.l d 19951) have been extens ively tested to explain 
the population of nearby white dw arfs ( Bergeron et al.l l2001 : 
Bergeron 2001; Be rgeron & Leggetn l2002) and use a very de- 
tailed description of the white dwarf atmosphere combined 
with the latest opacities (P. Bergeron 2005, priv. comm.). This 
is not the case for the models of Serenelli et al. and Hansen, and 
thus we should be careful in using their models quantitatively. 
Indeed, as can be seen from Fig. El their models do not repro- 
duce the observations of cool white dwarfs well. For instance, 
for the companion of PSR J0437-4715, which has a well- 
determin ed mass of 0.23 6 + 0.0 1 7 M an d distance of 1 39 ± 3 p c 
dvan Straten et al.ll200ll) . the models of ISerenelli et all d200ll) . 
while consistent with the observed B-R and R — I well, do not 
reproduce R — I and M R simultaneously. In contrast, the 0.2 M 
model of Hansen (2004, priv. comm.) does pass through the 
R — I, Mr point, but cannot reproduce both colours. It may 
be that both problems reflect uncertainties in the model atmo- 



spheres used by Hansen and ISerenelli et all d200ll) . It would 
be worthwhile to couple the evolutionary models of these au- 
thor s with the update d, very detailed atmospheric model of 
bergeron etldld 19951) . 

3.3. Cooling history and nature of the envelope 

Despite the uncertainty in the models and in the composition 
of the atmosphere, our observations show that the companion 
of PSR J075 1 + 1807 has cooled much more than expected if 
the amount of hydrogen was thick enough for significant resid- 
ual nuclear burning (Sect.[0. Indeed, the temperature is as ex- 
pected if no residual hydrogen burning occurred. F or instance, 
at the characteristic ag e of the pulsar, t = 7 .1 Gyr dNice et alJ 
2005), the 0. 196 M Q of lSerenelli etafl fcOOlh . which has a thin 
envelope, predicts a temperature of about 3200 K, which is 
roughly consistent with what is observed. With a pure helium 
atmosphere, a slightly colder temperature, of ~ 2500 K, is ex- 
pected, though t his is a less secure estimate due to uncertainties 
in the opacities (Hansen & Phinnev 1998a) 

The presence of a thin (or no) hydrogen envelope is not ex- 
pected, however, since thick envelopes are inferred for other 
optically identified companions in short-period systems (see 
Sect. [TJ. What could be wrong with this expectation? It was 
based on two theoretical ideas: (i) that below a certain criti- 
cal mass, no shell flashes occur and hydrogen layers will be 
thick; and (ii) that the companion mass monotonously increases 
with increasing orbital period. These assumptions appeared 
to be confirmed by the available data: for PSR J075 1 + 1807, 
with a perio d of 0.26 d. the c ompanion mass of 0.16-0.21 M 
(95% conf.: lNice etafll2005l) is similar to what is found for 
two other short-period systems with companions for which 
thick hydrogen envelopes are inferred, and less than the 
masses for longer period systems with thin-envelope com- 
panions. Specifically, PSR J 1 1 2+5307 (0.60 d, 0. 1 2-0.20 M ) 
and PSR J1909-3744 (1.53d, O.19-O.22M ) have thick en- 
velopes while PSR J0437-4715 (5.74 d, 0.20-0.27 M ) and 
PSR B1855+ 09 (12.33d, 0.24-0 29 Mp ) have thin envelopes 
(see Fig. |51and van Ker kwiik et all2005l and reference therein). 
Thus, while the uncertainties do not exclude that the compan- 
ion of PSR J075 1 + 1807 is so massive that it its envelope was 
diminished by shell flashes, the existing data make it unlikely. 

Two explanations for a thin envelope remain. First, there 
may be differenc es in metallicity among the progenitors of pul- 
sar companions. ISerenelli et alJ d2002l) studied the evolution 
of low-mass pulsar companions with sub-solar metallicity and 
found that, since the thermonuclear flashes are induced by the 
reactions of the CNO-cycle, the threshold mass between thin 
and thick hydrogen envelopes increases with decreasing metal- 
licity of the white dwarf progenitor. Thus, it may be that the 
companion of PSR J075 1 + 1 807 had a sufficiently higher metal- 
licity that it was above the threshold for shell flashes, while 
companions in other short-period systems had lower metallic- 
ity and hence were below the threshold, despite having higher 
masses. 

The next possibility is that the white dwarf was indeed 
formed with a thick envelope, which was subsequently re- 
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Fig. 3. The orbital period as a function of companion mass for 
a selection of low-mass binary pulsars outside globular clus- 
ters (either with P b < 5 d and M WDimin > 0.1 M©, or with P b < 
100 d and a secure companion mass determin ation). The data is 
comp i led from the A TN F Pulsar Catalogue jManchester et alJ 
20051 IStairsI ( 12004 and Ivan Kerkwiik et alJ J2005I) and refer- 
ences therein. Companion masses are either determined (dou- 
ble error bars, 95% confidence) or based on the assumption of 
a 1.35 M neutron star and an inclination of i = 90° (left er- 
ror bar, minimum mass), 60° (central symbol, median mass) or 
18° (right side, 5% probability that i is lower than this value and 
that the companion is heavier). Different central symbols indi- 
cate companions for which a thick or thin hydrogen envelope 
is inferred from optical measurements. The systems for which 
the neutron star mass is measured are indicated and connected 
by a gray line with their companions. The vertical grey lines in- 
dicat e the 1.35 + 0.04 M neutron star mass range determined 
bv lThorsett & ChakrabartvNl999l) . while the curved grey lines 
represent the theoretical relation bv iTauris & Savoriiia y999) 
between the white dwarf mass and the orbital period. 



moved by an action other t han shell flashes . Based on the upper 
limit o n the temperature of Lund gren et alJ lll996h . lErgma et all 
J200ll) already argued that the pulsar companion could not 
have the thick hydrogen envelope, and they proposed a scenario 
where part o f the envelope was removed by pulsar irradiation. 
lErgma et alJ found that irradiation driven mass-loss could re- 
move as much as 0.01 M from the thick hydrogen envelope 
(mostly while the companion is contracting following the ces- 
sation of mass transfer). 

A possible problem with the above suggestions, is that none 
predict the removal of the entire hydrogen envelope, while the 
observed colours seem most consistent with a pure helium or 
at least helium-dominated atmosphere. 



4. Irradiation by the pulsar? 

Above, we have treated the companion as if it were an isolated 
object rather than member of a binary system. Might the pres- 
ence of a relatively energetic pulsar influence our observations? 

The observed pulsar period and period derivative im- 
ply a spin-dow n luminosity Lsn = (2n) 2 IP/P 3 - 7.5 x 
10 33 hi ergs s" 1 (Eundgren e t alll995llNice et al.l2o"()5l) . where 
/ = 10 45 /45gcm 2 is the pulsar moment of inertia. For a 
2.1 M pulsar and a O.19M companion, the orbital separa- 
tion is a — 2.3 R , and, consequently, the irradiative flux of 
the pulsar wind incident on the companion is fi n = 2.1 x 
10 10 /45 ergs -1 cirT 2 . This is about twice the flux of the com- 
panion itself, / t h = o~T* s = 1.06 x 10 10 ergs -1 cm" 2 for 
Teff = 3700 K. Therefore, the presence of the pulsar and its 
irradiation may be important. 

Given the irradiation, one would expect the side of the com- 
panion facing the pulsar to be brighter than the side facing away 
from it. Thus, from Earth, the companion should appear faintest 
at phase 0.25 and brightest at phase 0.75 (using the conven- 
tion that at phase 0, the pulsar is at the ascending node). This 
is indeed seen in other pulsar binaries, with the black widow 
pulsar PSR B 1957+20 perhaps the most sp ectacular example 
(Ivan Paradiis et al!l988tlFruchter et al.ll988h . 

For star X, assuming a fraction 77 of the incident flux 
is absorbed and reradiated as optical flux, the flux from 
the bright side of the companion should be a factor 1 + 
f^/irr/Zth brighter (here, the factor | reflects projection effects). 
Observationall y, the inf erred values of 77 range from 0.1 to 0.6 
rtOrosz & van Kerkwijkfcoo'^ and references therein), and thus 
one expects a maximum change in bolometric flux by a factor 
1.13 to 1.8. For the 7?-band flux, the range is 1.2 to 2.2 (assum- 
ing it scales like a black-body spectrum, oc T b around 3700 K). 
We confirmed this u sing a de tailed light-curve synthesis model 
(described briefly in Stap pers et al.l 1999). 

For star X, n o effect is s een. U sing the PSR J075 1 + 1807 
ephemeris from iNice et ID (§005), we find that during the 
ESI /?-band observations the orbital phase ranged from 0.22 
to 0.25, while the 1996 LRIS /?-band images were taken at 
phases 0.86-0.90, and the 2005 LRIS images at phases 0.01- 
0.14 on January 7, and 0.77-0.93 on January 8. Thus, these 
observations span the orbital phases necessary to test for any 
modulation in brightness. I ndeed, using the inclination inferred 
from timing, i = 66^ deg INice et alJl2005l) . we find that dur- 
ing the ESI observations only 4 to 5% of the irradiated part 
of the companion surface was in view, while during the 1996 
LRIS observations is was 78% to 85%. As a consequence, 
we expect to see nearly the maximum change in brightness. 
Nevertheless, in Sect. 12.21 we found no significant variation, 
^lris - ^esi = 0.03 + 0.13; thus, to ~ 99% confidence, the 
variation is smaller than 0.3 mag, which implies 77 < 0.15. 

The lack of observed modulation could be taken to indicate 
that the irradiation is not very effective, e.g., because the albedo 
is large (i.e., 77 is small), the pulsar emission is non-isotropic, 
or the spin-down luminosity is overestimated. We believe these 
options are not ve ry likely (for a discussion in a slightly dif- 
ferent context, see Orosz & van KerkwiiklEool . which leads 
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us to consider the only alternative, that one of the assumptions 
underlying the above estimates is wrong. 

In particular, we assumed implicitly that irradiated flux is 
reprocessed and re-emitted instantaneously, i.e., transfer of flux 
inside and around the companion are assumed to have negligi- 
ble effect. For the companions of black-widow pulsars, this is 
reasonable, since for these relatively large objects, tides will 
have ensured synchronous rotation. Any flux transfer would 
thus have to be due to winds and/or convection, which plausi- 
bly happens on a timescale long compared to the thermal time 
of the layer in which the pulsar flux is reprocessed. 

The companion of PSR J075 1 + 1807, however, is well 
within its Roche-lobe, and tidal dissipation should be 
negligible. We can estimate its current rotation period 
from its prior evolution ? following the reasoning used by 
Ivan Kerkwiik & Kulkarni ( 1995) for the companion of PSR 
B0655+64. Briefly, during mass transfer, the companion filled 
its Roche-lobe and tides ensured the system was synchronised 
and circularised. Once mass transfer ceased and the companion 
started to contract to a white dwarf, however, the tides became 
inefficient, and the rotational evolution of the companion was 
determined by conservation of angular momentum. 

For our estimates, we split the total moment of inertia 
of the progenitor into two parts, one from the core, 7 core = 
^ore^core^core and one from the envelope, 7 env = kl m M em R^; 
here k is the radius of gyration and is the radius of the 
Roche lobe. After contraction of the envelope, one is left with 
a white dwarf with 7 WD = k^ /D M\ ND R^ ND . If we now assume 
that 7 core =i /wd and ignore differences in radius of gyration, 
conservation of angular momentum yields £2 rot /Q or b - 1 + 
M em R^/ MwdR\/ D - In reality, likely the envelope will be more 
centrally concentrated than the white dwarf, i.e., k em < &wd, 
and tidal dissipation will be important in the initial stages of the 
contraction. This will reduce the spin-up. On the other hand, the 
hot core of the progenitor will be larger than the white dwarf, 
i.e., /core > Iwd- In any case, it follows that unless the envelope 
mass is very small, the white dwarf should be significantly spun 
up. 

Model predictions for the envelope mass of helium-core 
white dwarfs differ. The O.196M model bv ISerenelli et alJ 
has an env elope mass of 6.7 x 10 3 M G (as given in 
Alth aus et all200i|). where a s a mo del of similar mass (Mwd = 
0.195 Mm) bv lDriebe etal] Jl998l) has one of 3.1 x 10~ 2 M G . 
Using these values, taking M W d = 0.19 Mo, Ry/D = 0.021 R G 
and /?l = 0.48 R Q , and ignoring differences in k, we find current 
rotation periods a factor 18-85 faster than the orbital period, or 
20 to 5 minutes. Given that thick envelopes seem inconsistent 
with the low observed temperature (Sect. [3j, the slower end of 
the range seems more likely. 

To estimate the timescale on which the pulsar flux is re- 
processed, we assume that the incident particles are predom- 
inantly highly energetic, and that they penetrate to, roughly, 
one Thompson optical depth. This corresponds to a column 
depth of N = 1.5 x 10 24 cm -2 , for which the thermal timescale 
t =t NkT/crT* s 1 min, where the numerical estimate is for 
T = T e g = 3700 K. This is shorter than the rotation peri- 
ods estimated above, suggesting that rotation may not be too 
important. On the other hand, our estimate is very rough. For 



instance, at one Thompson depth, the opacity at optical wave- 
lengths is much sma ller than unity for th e cool temperatures 
under consideration JSaumon et alJll994l) . Thus, the material 
likely radiates less efficiently than a black body, which would 
make the thermal timescale longer. Furthermore, the irradia- 
tion will change the temperature and ionisation structure of the 
atmosphere, further complicating matters. (Indeed, could this 
be the underlying cause for the fact that the colours deviate so 
strongly from those expected for a pure hydrogen atmosphere?) 
Finally, it might induce strong winds which equalise the tem- 
perature on both hemispheres (as is the case for Jupiter). 

5. Conclusions 

We have optically identified the white dwarf companion of the 
binary millisecond pulsar PSR J075 1 + 1807. We find that the 
companion has the reddest colours of all known millisecond 
pulsar companions and white dwarfs. These colours indicate 
that the companion has a very low (ultra-cool) temperature of 
r e ff ~ 3500 - 4300 K. Furthermore, the colours suggest that 
the white dwarf has a pure helium atmosphere, or a helium at- 
mosphere with some hydrogen mixed in, as invoked for the 
field white dwarf W D 0346 +246 which has similar colours 
(Qpp enheimer et all20oillBergerorj200ll) . 

Our observations are inconsistent with evolutionary mod- 
els, from which one would expect a pure hydrogen atmosphere. 
Indeed, as for other short-period systems, the hydrogen enve- 
lope is expected to be thick enough to sustain significant resid- 
ual hydrogen burning, leading to temperatures far in excess of 
those observed. It may be that the mass of the envelope was 
reduced due to she ll flashes or ir radiation by the pulsar, as was 
proposed bv lErgma et alJ J200ll) . 

However, we see no evidence for irradiation, despite the 
fact that the pulsar spin-down flux inpinging on the white dwarf 
is roughly double the observed thermal flux. Clues to what hap- 
pens might be found from more detailed studies of the spectral 
energy distribution, or more accurate phase-resolved photome- 
try. 

Finally, a deeper observation at infrared wavelengths would 
allow one to distinguish between the different atmosphere com- 
positions for the companion: for a pure helium atmosphere, 
black-body like colours are expected, while if any hydrogen 
is present, the infrared flux would be strongly depressed (as 
is seen for WD 0346+246). With adaptive optics instruments, 
such observations should be feasible. 
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